Binaries and core-ring structures in self- gravitating systems 
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Low energy states of self-gravitating systems with finite angular momentum are considered. A constraint 
is introduced to confine cores and other condensed objects within the system boundaries by gravity alone. 
This excludes previously observed astrophysically irrelevant asymmetric configurations with a single core. We 
show that for an intermediate range of a short-distance cutoff and small angular momentum, the equilibrium 
configuration is an asymmetric binary. For larger angular momentum or for a smaller range of the short distance 
cutoff, the equilibrium configuration consists of a central core and an equatorial ring. The mass of the ring 
varies between zero for vanishing rotation and the full system mass for the maximum angular momentum L ma x 
a localized gravitationally bound system can have. The value of L max scales as v/ln(l/a;o), where xo is a ratio 
of a short-distance cutoff range to the system size. An example of the soft gravitational potential is considered; 
the conclusions are shown to be valid for other forms of short-distance regularization. 

PACS numbers: 04.40.-b, 05.70.Fh, 98.10.+Z. 



I. INTRODUCTION 

Despite an abundance of rotating structures in the Uni- 
verse, the statistical mechanics of these systems is not com- 
pletely understood. The non-rotating self-gravitating systems, 
as shown by the Mean Field (MF) analysis and recently con- 
firmed by direct computer simulations, exhibit two phases 
(see, for example, QyllHlil and references therein): A high 
energy "uniform" phase where density contrast is small, and 
a low energy phase consisting of a diluted halo and a dense 
core. The structure of the core is determined by a type of a 
short-range regularization, which can vary from the exclusion 
in phase or coordinate space to a small-distance truncation of 
the interaction potential. At the same time, the structure of ro- 
tating self-gravitating systems, particularly the form and num- 
ber of cores or other dense objects in the low energy states, 
remains a subject of discussion. Configurations such as sym- 
metric and asymmetric binaries and rings QSHUllj a single 
spheroidal core with or without an equatorial ring fiol [Till , 
and a single "asymmetric" core, sliding along the container 
wall 1 12], have been found in rotating self-gravitating sys- 
tems. Among possible reasons for this plethora of equilibrium 
configurations the boundary conditions play definitely an im- 
portant role. As in the non-rotating case (see, for example, 
Ql), a confining box with reflecting walls (to conserve the 
energy and mass) is essential for the existence of any equilib- 
rium state. For a rotating system, it is natural to require an- 
gular momentum conservation as well and chose an axisym- 
metric boundary; usually a spherical container is considered 
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1 311 - However, the reflective boundary conditions 
seem to be the main reason for the apparent discrepancy be- 
tween the particle simulations and the MF results: While a 
true particle system inevitably evolves towards an asymmet- 
ric single-core configuration 1 12], the states with an arbitrary 
number of cores and even rings are predicted by the MF anal- 
ysis aS EzHH [01 • To obtain two- and multi-core states, the 
center of mass of the system is fixed in the center of the con- 
tainer |5[ |a LZl |8[ |9[ • While in the MF analysis the center of 
mass constraint is implemented by removing the dipole terms 
from the multipole expansions of density and potential, it is 



not clear how to fix the center of mass in a particle simulation. 
The total momentum and the center of mass position, unlike 
the angular momentum, are not conserved even in a spherical 
system, when the boundaries are reflecting. Due to the ergod- 
icity of a 3-dimensional Coulomb system, even a specially 
prepared, highly symmetric initial state (with a a centrally- 
symmetric coordinates and momenta distribution) will evolve 
towards the usually non-symmetric highest entropy configu- 
ration. For a rotating system, the most probable configura- 
tion consists of the single core sliding along the container wall 
a[3l- The reason why the state with a single core has the 
highest entropy (or the lowest free energy) is the following: 
When two or more initially separated cores merge, the gravi- 
tational potential energy decreases, which, for the fixed total 
energy, leads to a gain in the translational entropy of the halo. 
Intuitively, as there is no naturally occurring analog of a con- 
tainer wall on which the core may slide, this asymmetric state 
with a single core looks highly unphysical. 

If particle simulations seem unable to reproduce the con- 
figurations obtained by the MF methods and exhibit only a 
physically-irrelevant state with a single core, is there any other 
way to validate the MF results? The criteria of physical rele- 
vance, which the asymmetric state with a single core does not 
satisfy, are rather intuitive and can be formulated as follows: 
The physically relevant equilibrium states of rotating systems 
must be affected by the boundary conditions in the least possi- 
ble way. This minimal boundary condition effect is attained in 
the case of the core-halo states of non-rotating self-gravitating 
systems: If a container surrounding such system is removed, 
the halo will start to evaporate, while the core will remain al- 
most intact for a considerable time (see 1 14] for an estimate of 
core-halo thermalization rate). On the contrary, the rotating 
asymmetric state with a single core would undergo signifi- 
cantly more dramatic changes if the container were removed: 
The core, being no longer supported by the container wall, 
will escape ballistically with nothing left within the system 
boundaries. Similarly to the non-rotating case, the minimal 
boundary condition effect can be attained in rotating systems 
if cores (or other condensed objects such as disks or rings) are 
confined within the system boundary by gravity alone. In the 
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limit of the ground-state energy (or zero temperature), such a 
state will consist only of the gravitationally bound condensed 
parts, unaffected by the removal of the container at all. 

In this paper we consider such weakly interacting with 
container, or "physically relevant" states of rotating self- 
gravitating systems. To satisfy the relevance criteria suggested 
above, one needs to find core orbits in the presence of a halo. 
To simplify this task, we consider the limit of ground-state en- 
ergy (or zero temperature) in which the gaseous halo is con- 
densed into the cores, and there is no internal motion of core 
particles. The only motion in this limit is the macroscopic 
movement of cores, specified by the angular momentum con- 
straint. For sufficiently low energies, the state with the low- 
est ground-state energy is the thermodynamically equilibrium 
one (while the other mechanically stable states are thermody- 
namically metastable). Thus, to determine a core structure of 
the physically relevant low-energy equilibrium state, it is suf- 
ficient to find a gravitationally bound core configuration with 
the lowest energy of a given mass, size, and angular momen- 
tum. 

The number of candidates for the lowest energy core state 
can be significantly reduced using the following heuristic ar- 
gument. The energy of an ensemble of rotating cores is min- 
imized when the mass is concentrated into the largest core. 
This is so because for sufficiently small short-range cutoff, the 
total energy of a rotating self-gravitating system is dominated 
by the negative gravitational self-energy of the cores. The ab- 
solute value of gravitational self-energy of a core grows with 
the core mass faster than linearly for all reasonable forms of 
short-range cutoff. Hence, the energy is minimized by a con- 
figuration which consists of the principal core of the maxi- 
mum possible mass, while the remaining mass is distributed 
to carry the given angular momentum in the most energy- 
efficient way. For this reason the particle system in simu- 
lations always evolves towards the asymmetric configuration 
consisting of the single core which carries all the angular mo- 
mentum itself 1 12]. Without a constraint that the system has 
to occupy a finite volume, the lowest-energy state would have 
consisted of a core containing all but one particle with that 
particle having an orbit radius defined by the angular momen- 
tum constraint. Likewise, with the spatial localization con- 
straint, the most mass- and energy-efficient way to carry the 
given angular momentum is to put the smallest possible mass 
on a circular orbit of the maximum allowed radius. Being 
evident for a two-body system (see, for example, iTiUn . the 
energy efficiency of circular orbits can be seen from the fol- 
lowing argument: If a circular orbit is perturbed by adding 
a radial component to the velocity, the angular momentum is 
unchanged while the energy increases. Below we consider 
two possible configurations consisting of the principal core 
and the remaining mass on a single circular orbit: A central 
core with a ring of N orbiting cores and a binary, generally 
asymmetric. A symmetric binary is a limiting case of both 
families. Examples are abundant in the Universe and have 
been observed in the MF analysis 01 0] as possible equilib- 
rium or metastable configurations. Other mechanically sta- 
ble rotating core configurations that do not belong to these 
two families apparently have less mass concentrated into the 



largest core, and thus have higher energy. 

Using simple mechanics, we will derive that the choice of 
the lowest energy state depends on the range of a short-range 
(or high-density) regularization and the angular momentum. 
For an intermediate range of the small-distance regularization 
and small angular momentum, the binary state has the lowest 
ground-state energy, which confirms the results of |8]. In a 
limit of the vanishing range of the cutoff, or for the higher an- 
gular momentum, the core-ring state becomes the equilibrium 
one. The paper is organized as follows: after this introduction 
we define the model more formally. In Section III we com- 
pare the ground-state energies of two families of systems: A 
central core with a multicore ring, and an asymmetric binary. 
A Discussion and Conclusion section completes the paper. 

II. DEFINITION OF THE MODEL 

Let us now formally define the model. We search for the 
lowest-energy state of a system of M ^> 1 self-gravitating 
unit mass particles with fixed total angular momentum L lo- 
calized within a sphere of radius R. To make the problem 
analytically tractable, we will limit our consideration to the 
case when the range of small-distance (high-density) regular- 
ization is short. Hence the volume of all condensed objects 
in the ground state is considered to be negligible compared 
to the volume of the system. In addition, to make a clear 
distinction from the non-rotating case, we focus our atten- 
tion on the sufficiently high values of angular momentum L 
to exclude the configuration with the single spinning central 
core. Most of the analysis below is for a system of classical 
particles interacting via the attractive soft Coulomb potential 
— (r 2 + ^q) -1 / 2 (the gravitational constant is set to be unity). 
This simple form of short-range regularization is qualitatively 
equivalent to other forms of "softening", such as truncation of 
the Fourier or spherical harmonic expansions. For an inter- 
particle distance r smaller than the respective softening radius 
(given by r$, or by a wavelength of the highest untruncated 
harmonic function) all soft potentials tend to a harmonic os- 
cillator potential. Consequently, for all soft potentials a con- 
densed core with no particle motion (at zero temperature) is 
a point-like object. This is different from, for example, the 
ground state of a system of fermions or hard core particles 
with the finite core volume. Qualitative arguments will be 
given to show that conclusions made for the system with soft 
potential also hold for systems with other forms of short-range 
regularization. 

In addition to the notations r, E, and L for distance, energy, 
and absolute value of angular momentum, in the following 
we will also use the rescaled (universal) units for distance x, 
energy e, and angular momentum A (see, for example, (jj |3l 
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III. CORE-RING AND ASYMMETRIC BINARIES 



First consider a ring, i.e. a system of N identical cores 
(point masses in the case of a soft potential) of mass M/N 
each moving on a circular orbit of radius r < R. Assume 
that the softening radius tq is much smaller than the distance 
between the neighboring cores, and the interaction between 
cores can be considered as bare gravitational one. 

Mechanical equilibrium requires the centripetal accelera- 
tion of each core to be equal to the total force exerted on it by 
the others, 




FIG. 1 : A sketch of a central core - ring system. 
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where for large N the sum is replaced by an integral. This 
asymptotic expression turns out to be quite robust (less than 
4% off, even for N = 2), and the small numerical value of 
the constant C « 0.126 makes C negligible even for few-core 
systems. 

It follows from J2I3> that a binary system (JV = 2) local- 
ized within a sphere of radius R cannot have rescaled angular 
momentum higher than 



Xun = 1/V8« 0.35. 



(4) 



In a ring, the number of ring cores grows as N ~ exp(27rA 2 ) 
for large A. The increase of N, however, cannot continue in- 
definitely: For the interaction between cores to be similar to 
the bare gravitational potential — 1/r, the distance between 
cores should be larger than the softening radius ro. Hence, 
when the number of cores exceeds the corresponding limit 
N > N max 2irr/ro, the cores may still get closer to each 
other, while the maximum force acting on each core saturates 
at the corresponding to N max value. Consequently, the max- 
imum angular momentum of a self-gravitating system with 
radius not greater than R is 
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Qualitatively similar estimate for the maximum angular mo- 
mentum X m ax exists for systems with finite-size cores, such 
as those formed by ensembles of fermions or hard-core parti- 
cles. For such systems, the maximum number of cores N max 
corresponds to a merging of cores into a continuous (finite- 
volume) ring or torus. The radius of the body of the torus r c 
(which is of the order of the radii of cores before merging), 
serves the role of a cutoff parameter ro in (|5j- This conclu- 
sion can also be reached using the following argument: One 
can split a continuous ring into N max segments of size ~ r c 
and consider the interaction between them in the multipole 



expansion. The monopole-monopole interaction gives rise to 
the leading logarithmic term in l|3] |5}, while the higher-order 
terms produce only 0(r c /R) corrections. 

Let us now consider a core-ring system consisting of a sin- 
gle central core of mass M(l — a) and a ring of mass Ma 
consisting of N > 2 cores. This structure, sketched in Fig. [2 
resembles the planet Saturn with its ring. 

Similarly to the Eq. Q, the equation of motion for an or- 
biting core yields for the angular momentum of the system: 
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The total energy of the core-ring system consists of the grav- 
itational self-energies of the central core and orbiting cores, 
and the energy of macroscopic rotation: 



(1-a) 2 a 2 a 3 [a(lnA7(27r)-l] + l] 2 
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To find the ground state energy, this expressions must be min- 
imized with respect to a and N, 2 < N < 1/xq, taking into 
account the r < R constraint which has the form 
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For sufficiently small xq, the first term in @ is the dominant 
one. Hence, the minimum energy is reached when N is in- 
creased to its saturation value ~ l/%o to minimize the 
relative ring mass a within the range allowed by (|8jl. This is 
illustrated by the example presented in Fig. [2] The ground 
state energy of the core-ring configuration with the angular 
momentum A < IhN/2tt is 
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where a(A) is the minimal ring mass allowed by ^ and N = 
1/xq. 

The mass distribution of the ground states of the core- 
ring systems with other forms of cutoff is similar. The self- 
energy of a spherical core of m particles can be expressed as 
Ef e f f = -C HC m b / 3 and E[ elf = -C F m 7/3 for systems of 
hard-core particles and fermions, respectively (see 1 16] for the 
energy of self-gravitating fermion ball). The constants Cjjc 
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FIG. 2: Contour plot of e vs a and N, defined by the Eq. {7) and 
constraint JSJ for A = 0.5 and xq = 1CP 3 . The energy decreases 
(increases by absolute value) from dark to light, and reaches its min- 
imum at the largest possible N, N — 1/xo = 1000 and the smallest 
a allowed by {8j, a ~ 0.49. The black area corresponds to values of 
a and N which do not satisfy the localization constraint r < R 

and Cp, which depend on the hard-core radius and the num- 
ber of internal degrees of freedom, play the role of 1/Vrj in 
(|9}. For sufficiently large values of these constants (or equiv- 
alently, sufficiently small cores), the self-energy of the central 
core dominates over the self-energies of the ring cores and the 
energy of macroscopic motion. Hence, for a reasonable form 
of the short-range cutoff, the lowest energy core-ring configu- 
ration consists of the central core of the largest possible mass 
allowed by the localization constraint, and the orbiting ring of 
the system radius R which carries all the angular momentum. 

Now let us consider a binary, generally asymmetric |8|], as 
another candidate for the lowest energy configuration with a 
moderate angular momentum A < Abi„ (@}. Introducing an 
asymmetry parameter < S < 1/2 and the core separation 
r < 2R, denote the core masses and orbit radii as Mi = 8M, 
M-2 = (1 — 8)M, r% = (1 — S)r, and r2 = Sr, correspond- 
ingly. The total energy consists of the energy of rotation and 
the gravitational self-energy of the cores, and reads in rescaled 
units: 

2A 2 2x 

To find the ground-state energy, this expression has to be min- 
imized with respect to 5 while taking into account the con- 
straint that the radius of the largest orbit n should not exceed 
R, 

8 2 (l — 8) > A 2 . (11) 

For all reasonable values of xq (xq < 32A 2 /3) the gravita- 
tional self-energies of the cores dominates the total energy e 



in Eq. dlOi . Hence, the energy is minimized 

£bin = T {LI) 

2 2xq 

with the highest possible asymmetry, or the smallest possible 
value of 8 < 1/2 allowed by (II Q . The same is true for other 
forms of short-range regularization: Since the self-energy of 
a self-gravitating core grows faster than linearly with the core 
mass, a binary with the highest possible asymmetry will al- 
ways have the lowest energy. 

This confirms the conclusion made in 1 8] that for low en- 
ergy an asymmetric binary has a higher entropy than a sym- 
metric one, albeit with the restriction that orbits in a binary 
are gravitationally-supported only for A < \bi n ■ 

To compare the ground state energies of core-halo and 
asymmetric binary configurations, the cubic equations J8I1 li 
have to be solved for a given angular momentum A < A&j n , 
and the resulting a and 8 are to be substituted into the ex- 
pressions for energy ( 191 12i . For sufficiently small cutoff ra- 
dius xo -C 1, the self-energies of the cores, described by the 
first terms in ( 191 12b . give dominant contributions to the total 
energy. Qualitatively, when the cutoff radius is still not too 
small, so that |lna;o|/27r <C 1, the logarithmic term in (|8} 
is negligible and a « 8. In this case, due to the noticeable 
contribution of the self-energy of the smaller core, —8 2 /2xq, 
the binary system has the lowest energy. On the other hand, 
if | In xq I /2-7T ^S> 1, the logarithmic term dominates in (|8} and 
a < 8. In this case, the self-energy of the bigger central core 
of the core-ring system becomes smaller than the self-energies 
of both binary cores, and the core-ring system is the equilib- 
rium one. The A-dependence of the binary and core-halo ener- 
gies in the borderline case of xq w 10~ 7 is illustrated in Fig. [3] 
Hence for xq < 10~ 7 , the core-ring configuration is the equi- 
librium one for all values of angular momentum A, while for 
xq > 1CP 7 there exists a range of A < Am„ for which the 
binary system is the equilibrium one. Naturally, the energies 
of the core-ring and binary configuration coincide for A = 
where both configurations are reduced to just one central core. 

IV. DISCUSSION AND CONCLUSION 

In the previous sections we considered the low-energy equi- 
librium states of rotating self-gravitating systems and arrived 
at the following conclusions: 

To set a distinction between the physically relevant and ir- 
relevant states, we suggested that the core part of the relevant 
states has to be moving on gravitationally supported orbits 
and not come in contact with the system boundary. In the 
low-energy or zero-temperature limit, such equilibrium states 
become localized core-only states and can exist without a con- 
tainer. 

We found that, depending on the range of the short- 
distance regularization and the angular momentum, two pos- 
sibilities exist for the equilibrium configurations of rotating 
self-gravitating systems: For an intermediate range of short- 
distance cutoff and small angular momentum L < -J M 3 R/8, 
an asymmetric two-core binary configuration has the lowest 
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FIG. 3: Plot of the rescaled ground state energies of core-ring e c _ r 
(solid line) and asymmetric binary ej (dashed line) configurations vs 
angular momentum A for xo = 10 -7 . 



ground-state energy and thus is the equilibrium one |8|. For 
L > d M 3 R/8, or for a very small range of the cutoff and 
arbitrary angular momentum, the equilibrium configuration 
consists of a central core and an equatorial ring. 

The precise value of the cutoff range at which these two 
configurations have the same ground state energy depends on 
the nature of the cutoff. For the soft-core interaction potential, 
the crossover between two ground states takes place when the 
softening radius ro satisfies R/tq ~ 10 7 . For other forms 
of regularization, the role of ro is played by the radius of the 
body of the ring, or the radii of multiple "cores" that constitute 
the ring. The maximum angular momentum of a localized 
self-gravitating system scales as 



M 3 Rln(R/r ) 
2^ ' 



In this limit the central core vanishes and the system consist 
only of a ring. 

It would be desirable to check these conclusions by ei- 
ther particle simulations or the Mean Field analysis. Unfortu- 
nately, at this stage, none of these seems feasible. As shown in 
11211 and discussed in the Introduction, any finite-temperature 
particle simulations will lead to the physically irrelevant con- 
figuration with the single core sliding along the container wall. 
Another obstacle lies in the size of such a computation: To be 
able see the crossover between the binary and core-ring sys- 



tem one needs a ring consisting of « 10 7 cores, each consist- 
ing of at least one particle. A similar requirement on spatial 
resolution makes the Mean Field calculations very challeng- 
ing too. Indeed, even the existing calculations with fairly large 
cutoff |5j|g,|7|,|8|], may not have enough mesh points to resolve 
the core structure 11711 - However, the predicted crossover ra- 
tio -R/Vo w 10 7 is not at all astrophysically irrelevant and 
may be encountered in the planetary disks and even Saturn 
rings. Given that the typical radii of rings around the big plan- 
ets of the Solar system is of the order of 10 5 km 1181 . the 
dominant presence of ring particles of the size of 1 cm and 
less will make the core-ring configuration thermodynamically 
more stable than the binary one. 

The present consideration of the core-ring structure is rather 
schematic and at the current level does not allow us to ex- 
plain the fine structure of rings as gaps and spikes. Neither 
the Roche limit, which may set another bound on the exis- 
tence of the low-orbit self-gravitating binaries due to the tidal 
interaction, is considered. Yet even the present level of model- 
ing of core-ring structures permits to make important thermo- 
dynamical predictions, and a distinction between the existing 
core-ring models. For example, the central core-ring ground 
state configuration considered here may look similar to the 
core-ring configuration found for the system of fermions with 
fixed angular velocity in 11 ill . However, this similarity is only 
superficial: The ring observed in II ill is formed by particles 
which could not be supported gravitationally at the equator of 
the central core and were shedded off to the container wall, 
while in our case the ring is supported only gravitationally. 

Finally, it is interesting to speculate on up to what energy or 
temperature, the correspondence between the equilibrium and 
lowest ground-state energy configurations holds? Indeed, as 
the energy increases and the mass of the halo becomes com- 
parable or exceeds the combined mass of condensed objects 
(cores and rings), most of the angular momentum is carried 
by the halo. This may lead to a reduction and complete dis- 
appearance of a ring. This scenario looks especially plausible 
if the angular momentum is noticeably smaller than L max so 
a gravitationally-supported structure with smaller moment of 
inertia than that of an equatorial ring can carry it. In addition 
to numerical methods, this scenario can be analyzed analyt- 
ically by approximating a halo as a uniform gas, as done by 
Chavanis 1 19] in the case of a non-rotating state. We leave this 
analysis for future research. 
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